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Abstract. New exactly solvable rationally-extended radial oscillator and Scarf I potentials 
arc generated by using a constructive supersymmetric quantum mechanical method based 
on a reparametrization of the corresponding conventional superpotential and on the addition 
of an extra rational contribution expressed in terms of some polynomial g. The cases where 
g is linear or quadratic are considered. In the former, the extended potentials are strictly 
isospectral to the conventional ones with reparametrized couplings and are shape invariant. 
In the latter, there appears a variety of extended potentials, some with the same charac- 
teristics as the previous ones and others with an extra bound state below the conventional 
potential spectrum. Furthermore, the wavefunctions of the extended potentials are con- 
structed. In the linear case, they contain -I- l)th-degree polynomials with v ~ 0,1,2, .. 
which are shown to be Xi-Laguerre or Xi-Jacobi exceptional orthogonal polynomials. In the 
quadratic case, several extensions of these polynomials appear. Among them, two different 
kinds of {i/ + 2)th-degree Laguerre-type polynomials and a single one of {i^ + 2)th-degree 
Jacobi-typc polynomials with v = 0,1,2, .. . are identified. They are candidates for the still 
unknown Ar2-Lagucrre and A'2-Jacobi exceptional orthogonal polynomials, respectively. 

Key words: Schrodinger equation; exactly solvable potentials; supersymmetry; orthogonal 
polynomials 
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1 Introduction 

Since the pioneering work of Bargmann [1], there has been a continuing interest in generating 
rational potentials for which the Schrodinger equation is exactly or quasi-exactly solvable. In 
this respect, supersymmetric quantum mechanics (SUSYQM) [21 [3l HI [5l [6] , as well as the related 
intertwining operator method [7] and the Darboux algorithm [8], have attracted considerable 
attention (see, e.g., [9l[ini[IIl[I21[l3l[ia[ia[l6]). 

SUSYQM indeed provides us with a powerful technique for designing Hamiltonians with 
a prescribed energy spectrum from some known exactly solvable one. For one-dimensional 
systems, to which we are going to restrict ourselves here, the simplest approach based upon first- 
order intertwining operators uses a single solution of the initial Schrodinger equation, defining 
the so-called superpotential. The latter then furnishes a SUSY partner Hamiltonian, whose 
spectrum may differ by at most one level from that of the initial Hamiltonian. More flexibility 
may be achieved by resorting to nth-order intertwining operators with n > 1, constructed from n 



*This paper is a contribution to the Proceedings of the Eighth International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 21-27, 2009, Kyiv, Ukraine). The full collection is available at 
|http: / /www. emis.de/journals / SIGMA /symmetry2009.html | 



2 



C. Quesne 



solutions of the initial Hamiltonian [TTl HHJ [HI [201 [211 HSl [23] . These higher-order operators 
may turn out to be reducible into a product of lower-order ones or to be irreducible. 

In all cases, not only the SUSY partner potential and spectrum are explicitly known, but also 
the corresponding wavefunctions are given by some analytic formulae. The resulting expressions 
for the potential and the eigenfunctions, however, often look a good deal more complicated than 
those for the initial Hamiltonian (see, e.g., [M])- In such a context, to get a SUSY partner 
potential that is some rational function is not a simple task in general. 

An efficient method for getting potentials and wavefunctions expressed in terms of some 
elementary functions was proposed a few years ago and called algebraic Darboux transforma- 
tion |14l I15j . This procedure, whose idea goes back to an older paper [22], consists in using 
the solutions of the hyper geometric or confluent hypergeometric function whose log derivative 
is a rational function and which have been classified by mathematicians |26| . 

In a recent work [27], an alternative scheme has been devised. Known exactly solvable 
potentials are rationally extended by effecting a reparametrization of the starting Hamiltonian 
and redefining the existing superpotential in terms of modified couplings, while allowing for the 
presence of some rational function in it. The newly introduced parameters can be adjusted in 
such a way that the initial exactly solvable potential with reparametrized couplings becomes 
isospectral to a superposition of the same (but with its parameters left undisturbed) and some 
additional rational terms. Such a first-order SUSY transformation has been shown to be part of 
a reducible second-order one, which explains the isospectrality of the conventional and rationally- 
extended potentials with the same parameters. 

The purpose of the present paper is twofold: first to demonstrate the efficiency of this 
method in producing rather complicated exactly solvable rational potentials, whose existence 
would otherwise have remained unknown, and second to examine in some detail the polynomials 
appearing in the rational potential bound-state wavefunctions. Among these polynomials, the 
simplest ones belong to two classes of exceptional orthogonal polynomials, which have been the 
topic of some recent mathematical study |28[ I29j . We plan to show here that our SUSYQM 
approach provides us with a convenient way of constructing and generalizing such polynomials. 

This paper is organized as follows. In Section [21 the method is first applied to rationally 
extend the radial oscillator potential and to study the associated Laguerre-type polynomials. 
The generalization to the Scarf I potential and to the corresponding Jacobi-type polynomials is 
then dealt with in Section [3l Finally, Section [J] contains some final comments. 

2 Radial oscillator potential 

The radial oscillator potential 



where to and I denote the oscillator frequency and the angular momentum quantum number, 
respectively, is defined on the half- line < x < oo. It supports an infinite number of bound 
states, whose energies are given by 



The associated wavefunctions, vanishing at the origin and decaying exponentially at infinity, 
can be expressed in terms of Laguerre polynomials as [3l [30] 










'U>X' 



,2 



(2.3) 
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where 



wxia+f) / 2 1/! 



1/2 



Before applying SUSYQM to such a potential, let us briefly review some well-known proper- 
ties of this method [H El H El E] . 

In the minimal version of SUSY, one considers a pair of first-order differential operators 

A = -^ + W{x), = -A. + w{x), (2.4) 

and a pair of factorized Hamiltonians 

=iti = — ^ + -S, =iit = --^ + F(-)(x)--E, (2.5) 

dx'^ dx^ 

where E is the so-called factorization energy. To E one can associate a factorization func- 
tion (j){x), such that H^~^^(f)(x) = E(f>{x). The superpotential is given in terms of (t>{x) by the 
relation W{x) = —<j)'{x)/<j){x) and the two partner potentials V^^\x) are connected with W{x) 
through the equation 

V'^^\x) = W^{x)tW'{x)+E, (2.6) 

where a prime denotes derivative with respect to x. 

Whenever E and (j){x) are the ground-state energy E^^ and wavefunction ipl^\x) of the 

initial potential y'^+)(x), i.e., H^^^(l){x) = H^^^ipl^\x) = 0, the partner potential V^~'){x) 
has the same bound-state spectrum as V^~^\x), except for the ground-state energy which is 
removed (case i). For E < E^\ in which case (j){x) is a nonnormalizable eigenfunction of H^^\ 
the partner V^~\x) has the same spectrum as V^~^\x) if (j)~^{x) is also nonnormalizable (case ii) 
or has an extra bound-state energy E below E^'^\ corresponding to the wavefunction (j)~^{x), if 
the latter is normalizable (case in). 

The standard SUSY approach to the radial oscillator potential corresponds to case i. Deletion 
of the ground-state energy Eq^ yields a partner potential Vi^i{x) of the same kind, this meaning 
that the radial oscillator potential is shape invariant [31]. In this case, 4>{x) = x''^^ ex.p{—^uix'^) 
or, in other words, 

W{x) = -ujx - (2.7) 
2 X 

We now plan to modify such a superpotential in order to generate rationally-extended radial 
oscillator potentials. 

2.1 Extended classes of rational radial oscillator potentials 

Motivated by the form of the standard superpotential ()2.7I) and by the fact that z = is 
the basic variable appearing in the Laguerre polynomials that control the bound-state wave- 
functions (j2.3p . let us make the ansatz 

N b dq/dx 

W{x)=ax+ ^ — , (2.8) 

X g 

where a, h are two constants and g is some polynomial in z (or simply x^). Our purpose is 
to show that it is possible to choose the constants a and 6, as well as those appearing in g, in 
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such a way that V^^\x), as defined in (j2.6p . continues to belong to the radial oscillator family 
up to some change in I, while its partner V^~\x) differs from Vi{x) in the presence of some 
additional rational terms, i.e., 

= V,{x), = Vi,ext(x) =Vl{x)+ Vl,rM- 

To be physically acceptable, such additional terms Vi^j.a,t{x) must be singularity free on the 
half- line < j; < oo. 

In the following, we carry out this programme for linear and quadratic polynomials g{x'^). 

2.1.1 Linear case 

Let us assume that 

g{x') = + c, ^ = 2x, 

where c is a constant satisfying the condition c > in order to ensure the absence of poles on 
the half-line. 

On inserting the ansatz (j2.8p in (j2.6p . we obtain 



y = a^x^ + "-^^ + 2^^^^^^^ + a{2h -5) + E. 



Elimination of the term proportional to (x^ + c) ^ and of the additive constant from F (x) 
can be achieved by choosing 

c=^^^, E = -a{2b-5). 
The partner potential is then given by 



x"^ ~^ x^ + c (x^ + c)2 ~'~ 
and its first two terms coincide with Vi{x), defined in (12. ip . if the two conditions 

a2 = la;2, b{b - I) = l{l + 1) 

are fulfilled. These two restrictions lead to four possible solutions for {a,b): {a,b) = {^,1 + l), 

(f,-0, + (-1,-0- 

As a result, c becomes c = ±(2/ -I- l)/uj, where only the upper sign is allowed. Since this 

is achieved either for {a,b) = + or for {a,b) = we arrive at two distinct 
solution^, 

yW(x) = v^^^ix), y(-)(x) = mx) + ^ - , '^'["111,, ± (2.9) 

wx^ + 2/ -|- 1 (wx^ + 2/ + 1)^ 

corresponding to 

Wix) = ±-u;x ± i±i±i - -^^^ , (2.10) 

E = -io{l + ^T2) , (2.11) 

</>(x) = x^('+5)-5(u;x2 + 21 + l)e=f3'^^', (2.12) 



^Strictly speaking, such a picture is only valid for positive I values, since for / = we are only left with the 
first solution. 
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where we take either all upper or all lower signs, respectively. We can therefore go from a con- 
ventional radial oscillator potential to an extended potential Vi^cxti^) in first-order SUSYQM 
provided we start either from Vi^i{x) or from Vi-i{x) and, apart from some irrelevant additive 
constant, the result of the extension is the same in both cases. 

Since the factorization energy E, given in (j2.1ip . is smaller than the ground-state energy 
of V^+\x), namely sj^^^ = cj (/ + | ± l), it is clear that we must be in case ii or case in of 
SUSYQM. On examining next the behaviour of the factorization function, it becomes apparent 
that it is the former case that applies here, namely V^~^^ (x) and V^~^ (x) are strictly isospectral. 
We indeed observe that for the first choice corresponding to upper signs in (j2.12|) . cj) decays 
exponentially for x oo, but grows as x"'"^ for x — > 0, while for the second choice associated 
with lower signs, it vanishes for x — > 0, but grows exponentially for x — > oo. 

We conclude that the addition of the two rational terms Auj{ujx^ + 21 + l)-i and -8uj{2l + 
l){iox'^ + 21 + 1)"^ to V/(x) has not changed its spectrum, which remains given by equation (|2.2p . 



2.1.2 Quadratic case 
For 

g{x^) = x^ + cx^ + d, M^!)=2x(2x2 + c), 

dx 

where c and d are two constants, the picture gets more involved because the singularity-free 
character of the extended potential on the half-line can be ensured in two different ways: the 
quadratic polynomial g{x'^) may have either two real and negative roots or two complex conjugate 
ones. In association with these two possibilities, we will have to distinguish between the two 
cases 

c> 0, 0<d<—, (2.13) 



and 

.2 

4 

On proceeding as in the linear case, it is straightforward to get 



d>^. (2.14) 



l^W (x) = a^x^ + ^ + ^2iac-2b + 3y + 4ad-i2b-l)c _ 



X 

leading to the conditions 



26-3 , (26-l)(26-3) ^ . , ^ 

, d = ^ ^, E = -a{2b-9). 

a 4a^ 



Furthermore 



2 2 ^(&-l) 4(2x^-c) 8{c^-id)x 

V'^ >{x) = a^x^ + ^ „ ^ + + „ ^ + 2a 

x* + cx^ + d [x^ + cx^ + dy 

yields the same four possibilities for (a, 6) as before. To go further, we have to distinguish 
between (l2l3D and (ITlil) . 

In the real case, the condition c > can be achieved by imposing a = ^uj, b = I + 1 with 
/ > 0, or a = —^oj, 6 = 1, or else a = ~^oj, b = —I. The next restriction d > 0, being equivalent 
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to (26 — l)/a > 0, discards the choice a = 6 = 1. FinaUy, d < ^ turns out to impose 

a < 0. Hence we are left with a single solution corresponding to 



--UJ, o = -t, c= — , d = 5 — , 7 = 2/ + 3. 

2 a; 



As a result, we get for / > the pair of partner potentials 
V<+>(i) = F,_i(i), 

Vi-\.) = V,i.) + , ""If- + ,^ ^ - 7.2/ + 3, (2.15) 

(u;x^ + 7)^ — 27 [(wx"^ + 7)^ — 27]"^ 

corresponding to 

1 / 4wx(w2;2+7) 

VV (x) = LOX 



2 X (wx^ + 7)2 — 27 ' 
^ = -cj(/ + f) , 

(/.(x) = x'[(wx2 + 7)2 - 27]e3^^'. (2.16) 

2 

Turning ourselves to the complex case, we observe that the single condition d > ^ is equiv- 
alent to 6 > |, which imposes 6 = Z + 1 with I > 0. The sign of a remaining undetermined, we 
have actually here two distinct solutions according to the choice made, namely 

1 27 7(7 + 2) 

a = ±-uj, h = l + l (/>0), c = ± — , (i=^Al_^^ ^ = 21-1. 

2 u> uj^ 

They lead to 

Vl+>(i) = V,+i(i), 



2 „A M'yuj'^x^ 



(wx^ ± 7)2 + 27 [(wx2 ± 7)2 + 27]2 

for / > 0. Correspondingly 

1 / + 1 4a;x(wx2±7) 
W{x) = ±-ujx + 



iba;, 7 = 2/-l, (2.17) 



2 X (u;x2 ± 7)2 + 27 ' 

E = ^u{l-l), 

(j){x) = x~^~\ujx^ ± 7)2 + 27]e^3^^'. (2.18) 

In comparison with the linear case which has provided us with a single rationally-extended 
potential, the quadratic one reveals itself much more flexible since three distinct rationally- 
extended potentials with / > have been obtained. It is also worth observing that the angular 
momentum of the starting potential has now a strong influence on the resulting extended one. 

Furthermore, we note that in all three cases, the factorization energy is smaller than the 
ground-state energy of V^~^\x)^ as in the linear case. As far as the factorization function is 
concerned, however, a distinction must be made between the function in (j2.16p or that with 
upper signs in (|2.18p . on one hand, and that with lower signs in the latter equation, on the 
other hand. For the first two functions, the situation is similar to the previous one leading to 
strictly isospectral partner potentials (case ii), whereas, for the third function, we observe that 
(/)~^(x) is normalizable on the half-line and is therefore the ground-state wavefunction of the 
corresponding V^~\x) potential with an energy eigenvalue -Eg ^ = Eq^^^ — 6u; = w (/ — |) . 

In the following, we shall refer to the potentials (|2.15|) and (|2.17|) with upper /lower signs as 
case I, II, and III potentials, respectively. 
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2.2 Determination of wavefunctions 

The factorized forms (|2.5p of H^^^ and H^~^ imply the existence of the intertwining relations 

El SI EE] 

which enable us to determine the wavefunctions of H^~^ from those of H^^^ with the same 
energy (or vice versa). Hence, on starting from H^^^ 'il)\j^\x) = £yip\^\x) with = ^ — E, 
we can write 

^Pi-Hx) = ^Aipi+\x), u = 0,1,2,..., (2.19) 



4 \x) ^ cl)-\x), 4J^{x) = ^A4+\x), 1^ = 0,1,2,..., (2.20) 



or 



according to whether case ii or case in of SUSYQM applies. 

We will now proceed to use equations (|2.19p and (j2.20p to determine the wavefunctions of 
the rationally-extended radial oscillator potentials constructed in Section [2.11 

2.2.1 Linear case 

In the linear case, the result for Tpi \x) should be the same whether we start from xpi^\x) = 

'ij:u~^^\x) or from ^pl^\x) = ipu Considering the latter case to start with, we can rewrite 
the operator A as 

/ / d 1 I 2 \ 1 2 

A = \2bJZ — ; , Z= —OJX . 

\dz 2 2z 2z + 2l + lJ' 2 

Here we have used equation (|2.4p with W{x) given in (j2.10p . where we assume the lower sign 
choice. On applying this operator on the right-hand side of equation ()2.3p with / — 1 substituted 
for /, and employing equation (j2.19p . we get 

^^.«;^^aS"'4"-"(.), (2.21) 



ujx^ + 21 + 

where A/'i'^'* = Nu~'^\ e,^ = 2u; (z^ + / + |) , a = Z -|- ^, and we have defined 

The action of the first-order differential operator o'f^^ on the Laguerre polynomial {z) 
can be inferred from known differential and recursion relations of the latter [32]. The result can 
be written as 

d'r'^L^r'\z) = L^:l{z), 

where L^'^^;^(z) is a (z^ + l)th-degree polynomial, defined by 

L^:l{z) = -{z + a + l)4")(z) + 4"_\(z) (2.22) 
in terms some Laguerre ones. 
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The partner potential wavefunctions (j2.2ip can therefore be expressed in terms of L^^^-^{z) as 
=^^~ \/+ll + i ^^^ti'^ (^--') (2.23) 

where 



^2'"i(. + / + i)r(. + / + i: 



Similarly, for the other choice ip\^\x) = ipy'^^'^ {x), we can write 

^ ^"V^^ a;x2 + 2Z + l 2 K l^J, 
where 7Vi+^ = 7Vi'+^\ e,. = 2u; (z^ + ^ + i) , a = Z + i, and 

of ^ = (z + a) (^^^ + a + l^-z. 
On taking into account that this time 

we arrive (up to some irrelevant overall sign) at the same wavefunctions (j2.23p as in the previous 
derivation, as it should be. 

At this stage, some observations are in order. In |28^ 129] . it has been shown that the family 
of {v + l)th-degree polynomials l'^^^-^{z), i/ = 0, 1, 2, . . ., defined in equation (|2.22p . arises in a 
natural way when extending Bochner's theorem on the relations between classical orthogonal 
polynomials and solutions of second-order eigenvalue equations [33] by dropping the assumption 
that the first element of the orthogonal polynomial sequence be a constant. It is actually 
one of the two complete orthogonal sets of polynomials with respect to some positive-definite 
measure that can be constructed by starting instead with some linear polynomial (the second 
one being considered in Section [3. 2p . For this reason, the exceptional polynomials i}i^-^{z) have 
been called Xi-Laguerre polynomials [281 129j . They satisfy a second-order eigenvalue equation 
with rational coefficients and are normalized in such a way that their highest-degree term 
is equal to {—lY~^^z^~^^/v\ (as compared with {—lYz'^/v\ for L^\z)). On combining their 
definition (j2.22p with the recursion relation of classical Laguerre polynomials, one can write 
them as linear combinations of three Laguerre polynomials with constant coefficients, 

L^+i(z) = {y + l)LZ^i{z) -2{v + a + l)LZ{z) + {v + a + (2.24) 

A first few polynomials are listed in Appendix |Al where they may be contrasted with the 
corresponding classical Laguerre ones. 

The extended radial oscillator potential (j2.9p and its corresponding wavefunctions (|2.23p were 
derived for the first time in [34] by using the point canonical transformation method [35] and 
the differential equation satisfied by L^'^^(z). In the same work, the potential was also shown 
to be shape invariant with a partner given by V/+i^ext(2;)- Here we have adopted another type 
of approach based on first-order SUSYQM, which has the advantage of defining the exceptional 
polynomials through the action of some first-order differential operator on a classical Laguerre 
polynomial. As we now plan to show, this procedure can be easily extended to the quadratic 
case for which it will lead to some novel results. 
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2.2.2 Quadratic case 

As proved in Section 12.1.21 the quadratic case leads to three distinct new rationally-extended 
potentials, of which the first two are similar to that obtained in the linear case (case ii of 
SUSYQM), while the third one is not (case Hi of SUSYQM). 

By proceeding as in the linear case and defining z = ^u>x^ and a = ^ + ^, we arrive at the 
following three first-order differential operators 

d["^ = + 2{a + l)z + a{a + 1)] (^-^ - 1^ - 2{z + a + 1), 
Oj"^ = [z^ + 2{a - l)z + a{a - 1)] + a + - 2z{z + a-l), 

of ^ = [z^ - 2{a - l)z + a{a - 1)] ( z-^ - z + a + l) - 2z{z - a + 1), 



dz 

leading to three distinct families of Laguerre-type polynomials, 

= -[z^ + 2{a + 2)z + {a + l)(a + 2)]4")(z) + 2{z + a + l)4-i(2)> (2-25) 

0('^)4"+i)(z) = (z. + a-l)Lg^2(^) 

= {{u + a- l)[z^ + 2az + {a - l)(a + 2)] + 2{a - l)}L^^\z) 

- 2(i/ + a)(z + a- l)4"\(z), (2.26) 
0(°)4"+^)(z) = (z. + 3)Lg+3(z) 

= {-z^ + {2u + 3a- 3)z^ - [2{2a - 3)iy + 3a{a - l)]z 
+ {a- l)[2(a - + a{a + 1)]}4"H^) 

- (z. + a)[z^ - 2{a - 2)z + {a - 1)(q - 2)]4-i(-z), (2-27) 

for cases I, II, and III, respectively. 

The corresponding rationally-extended potential wavefunctions can be written as 

=^-'' (.x^ + 2i-1)^' + 2(2i-l) ^t1i (i-^) e-i"'^ (-e II). 
*-(^' =^''- \..^-2lf^Y + 2i2l-l) tM iW) (-e III), 



where 



/ 1+3 \l/2 

AAi-) = - ^ (easel), 



AA(-) = ^^-^ (case II), 



3 \ 1/2 



Ml+i = [ r.y'^ ; i (case m), 

and v = 0,1,2, .. .. 
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Let us examine in more detail the Laguerre-type polynomials obtained in the first two cases. 
From physical considerations related to general properties of SUSYQM, it follows that the 
wavefunctions if^y \x), v = 0,1,2,..., form two complete, orthonormal sets in the Hilbert 
spaces associated with the corresponding rationally-extended potentials. This suggests that 
a similar property should be valid for L^^j^jl-^) and -£2"i+2(-^) with respect to some appropriate 
positive-definite measure. From the additional fact that the lowest-degree polynomials are 
quadratic in z (see Appendix lAj) . we infer that L^^l_^_2{z) and l'^1_^_2{z) are good candidates for 
the still unknown X2-Laguerre polynomials. Note also that by analogy with classical and Xi- 
Laguerre polynomials, in defining them in (|2.25p and (j2.26p . we have chosen their normalization 
in such a way that their highest-degree term is equal to (— l)'^+^z'^+^/z^!. Furthermore, it is 
straightforward to show that the counterparts of equation (j2.24p are now linear combinations of 
five Laguerre polynomials with constant coefficients, 

^g+2(^) = + 2)(z. + l)L%{z) - A{v + l){v + a + 2)<\(z) 

+ 2{v + a + 2)(3z^ + 2a + 2)4") (z) - 4(i/ + a + 2)(i/ + a)4-i(-2) 
+ {u + a + 2){v + a- l)L^^l2{z) 

and 

4^+2(^) = {v + 2){v + l)L%{z) -A{v + l){v + a + l)4^i(z) 

+ 2{v + a + l)(3z^ + 2a + 1)4"^ (z) - 4(i/ + a + l)(i/ + a)4-i(-2) 
+ (i/ + Q + l)(i/ + a)4-^2(-2)- 

Since for case III, the wavefunctions have to be supplemented with the function 

(f)~^{x) in order to form a complete set of orthonormal wavefunctions for the corresponding 
rationally-extended potential (see equation (j2.20p ). it is rather clear that the Laguerre-type 
polynomials l!^^'1_^^{z) will not form a complete orthogonal set of exceptional polynomials. One 
of their properties is however worth pointing out. Although from their definition ()2.27p one 
would expect that their expansion in terms of classical Laguerre polynomials would contain 
seven terms, one actually gets only five of them, i.e., 

4^+3^ = (^+ 2)(^+ l)L%{z) - A{v+ 2){v+ l)L%{z)+ 2{v+ 1)(3^+ a + 5)4^1 (z) 
- A{v + l){v + a + l)L^^\z) + {y + a + + a)4-i(^), 

the coefficients of L^^'^^gl-^) and l'^^}_,^{z) vanishing identically. Observe that in (j2.27p . the polyno- 
mials have been normalized in such a way that their highest-degree term is {—lY^'^z'^~^'^/[{v + 
3)z.!]. 

In cases I and II, the extended potential ground-state wavefunction ip\^ \x), expressed in 

terms of 2 ^ (i'^^^) o^" ^2 2 ^ (l"^-^^)' respectively, can be used to prove that the correspon- 
ding potential Vi^exb{x) is shape invariant, its partner being given by Vi+i^ext{x) (see AppendixlD]). 
This generalizes to the quadratic case a property already demonstrated for the linear one in |34j . 

3 Scarf I potential 

The Scarf I potential [3] 

Va,b{x) = [A{A - 1) +B'^]sec^x- B{2A - l)secxtanx. 
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defined on the finite interval —^<x<^, depends on two parameters A, B, such that < B < 
A — 1. It has an infinite number of bound-state energies 

Ei^^ = {A + uf, u = 0,1,2,..., (3.1) 

depending only on A, and with corresponding wavefunctions 

V;(^.^)(x)=AA(^'^)(l-sinx)^(^-^)(l + sinx)5(^+^)pi^-^-^'^+^^ z^sinx, 

.r{A,B) _ ( {2A + 2u)u\T{2A + u) \ ^'^ 

\2^^V{A-B + v + \)T{A + B + u + \)) ' ^^'^^ 

vanishing at both end points and expressed in terms of Jacobi polynomials. 

In the standard SUSY approach corresponding to case i, deletion of the ground-state en- 
ergy E^'^ yields another Scarf I potential Va+i,b{x) (hi other words the potential is shape 
invariant [31]). In such a case, 0(x) = (1 — sinx)('^~^-*/^(l + sinx)^"^^^-*/^ and 

W{x) = A tan x — B sec x. 

As in Section [21 we will now proceed to modify such a superpotential in order to generate 
rationally-extended Scarf I potentials. 

3.1 Extended classes of rational Scarf I potentials 

Let us assume a new superpotential of the form 

dg/dx 



W{x) = a tan x + 6 sec x 



9 



where a, b are two constants and g is some polynomial in z = sinx. We now demand that 
V^^^ (x) belongs to the Scarf I potential family up to some reparametrization of coefficients and 
that y(-)(x) differs fr om Va,_b(x) in the presence of some rational terms, i.e., 

y W(x) = Va',B'{x), V^-\x) = Va,bM^) = Va,b{x) + Kl,B,rat(x). 

Here VA,B,rat(x) should be a singularity- free function on the interval — ^ < x < ^. 

3.1.1 Linear case 

On assuming 

■ dg(sinx) 

g(smx) = smx + c, ; = cosx, c > 1, 

dx 

and proceeding as in Section 12.1.11 we arrive at the conditions 

E={a + lf, (3.3) 



2a + 1' 

and 



a{a + l) + b^ = A{A-1) + B^, {2a + l)b = -{2A - 1)B, (3.4) 

coming from and V^~\ respectively. 

By successively using equations (|3.4|) and (|3.3|) . we get {a,b) = (=F-B — i, ±A =F ^) or 
{^A - i ± ^,±B) and c = -{2A - 1)/{2B) or -2B/{2A - 1). Since only the former choice 
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for c is compatible with the condition |c| > 1, we are only left with (a, b) = {^B -^,±A^^), 
leading to 

V^+\x) = Va,b±i{x), 

.^ 2(2^-1) 2[(2A - 1)2 - 4^2] ^ ^ 

2A — 1 — 2B sm x [2A — 1 — 2B sm xj^ 

Correspondingly, 

TX./ ^ f A , f . A 2BCOSX 

= tB tanx ± M secx + 



2/ V 2/ 2yl- 1 - 2Bsin2;' 



= (1 + sinx)'F5(^+^--5±l)(l - sinx)^5(^--S-|=Fi)(2^ - l - 25sinx), 

where we take either all upper or all lower signs. 

It is straightforward to see that E and (/>(x) satisfy the conditions for strict isospectrality 
(case ii of SUSYQM). Since, on the other hand, the B value does not influence the Scarf I energy 
eigenvalues, it follows that the rationally-extended potential V^~\x) has the spectrum (13. ip . 
The results are therefore very similar to those of Section 12.1. It we have two possibilities for the 
starting reparametrized Scarf I potential, but a single one for the extended potential. 



3.1.2 Quadratic case 

With the assumption 

/ • N ■ 2 • 7 dg(smx) _ . , 

(7(sinxj = sm x + csmx + a, ; = cosx(2sinx + cj, 

ax 

the discussion gets rather involved. Firstly, one has to impose that the function g{z) has either 
two real zeros outside the defining interval — 1 < z < 1 or two complex conjugate ones and, for 
the former possibility, one has to distinguish between two zeros on the same side of the interval 
or one zero on the left and the other on the right. As a result, the parameters c and d may vary 
in three distinct ranges, 

l<\c\— l<d<— (case la), d < — |c| — 1 (case lb), 



and 

d > — (case 2). 



c2 



Secondly, one has to translate these conditions on c, d into restrictions on a, b, taking into 
account that the two sets of parameters are related to one another through the equations 

_ 46 ^_ 462 -(2a + 3) 



2a + 3' 2(a + l)(2a + 3)' 

coming from V^'^\ Thirdly, the resulting conditions on a, 6 have to be combined with the links 
(a, 6) = (=F-B - i, ±A T ^) or {=fA - ^ ± i, ±B) imposed by V^-l 

This leads to several choices for (a, 6) and for the domain of variation of A and B, which 
are summarized in Table [TJ Apart from the three options marked with an asterisk (and to be 
referred to as cases I, II, and III), most of the possibilities are characterized by a very limited 
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Table 1. Supcrpotcntial parameters for rationally-extended Scarf I potentials in the quadratic case. 



Case 


a 




6 


Conditions 


la* 


-B 


1 

2 




1< S < ^- 1 


lb 


-B 


1 

A 


Z 


|<A<2, }.<B<A-l 

Z ' Z 


lb 


-B 


_ 1 

2 




^ > 2, \ <B <l 

— ' Z 


2 


-B 


1 


z 


1 < A< 1, ^-A<B<A-l 

4 z ' z 


2 


-B 


1 

2 


2 


^>i, o<s<i 


2* 


B - 


1 

2 




< s < ^- 1 


2 


-A 




B 


1<A<|, 0<S<^-1 


2 


-A 




B 


|<A<|, 0<-B<-^+| 


2* 


-A 




B 


< < yl- 1 



range for at least one of the parameters A, B. For simplicity's sake, they will not be considered 
any further. 

For the selected cases, the two partner potentials can be written as 

^(+)(^) = Va'M^)^ y'^-H^) = VaA^) + §^ + ^> (3.6) 

where 

A' = A, B' = B + 1, 

iVi(x) = -A[{2A - 1){2B - 1){2B - 2) sinx + 2{2A - if - {2B - 2f{2B + 1)], 
N2{x) = -8{2B - 2){2A - 2B + 1){2A + 2B - 3) 

X [2(2^ - 1){2B - 1) sinx - {2A - if - 2B{2B - 2)], 
D{x) = {2B -1)[{2B -2)smx - {2A-l)f - {2A-2B + 1){2A + 2B -3) (3.7) 

for case I, 

A' = A, B' = B-1, 

Ni{x) = -A[{2A-1){2B + 1)(2B + 2)sinx - 2{2A- if - {2B + 2f{2B - 1)], 
N2{x) = 8{2B + 2){2A - 2B - 3){2A + 2B + 1) 

X [2{2A - 1){2B + 1) sinx - {2A - if - 2B{2B + 2)], 
D{x) = {2B + 1)[{2B + 2) sinx - {2A - l)f + {2A - 2B - 3){2A + 2B + 1) (3.8) 

for case II, and 

A' = A+1, B' = B, 

iVi(x) = -8[B{2A - 2){2A - 3) sinx - A{2A - 3f + AB^], 
N2{x) = 8{2A - 3){2A -2B- 3){2A + 2B-3) 

X [AB{2A - 2) sinx - 45^ - {2A - 1){2A - 3)], 
D{x) = {2A - 2)[{2A - 3) sinx - 2Bf + {2A - 2B - 3){2A + 2B-3) (3.9) 
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for case III. Correspondingly, 

E={B-^)\ (l){x) = (l-sinx)3(^-^-i)(l+sinx)-|(^+^)Z?(x), 
E = (5 + |)^ = (1 - sinx)-^(^"-^)(l + sinx)i(^+^-i)L>(x), 

and 

E = {A-2f, cj){x) = (l-sinx)-3(^-^)(l + sinx)-5(^+^)L)(x), 

from which it follows that cases I and II are characterized by strict isospectrality (case ii of 
SUSYQM), whereas in case III, (^~^(x), which turns out to be normalizable on (— f ) and to 
vanish at both end points, is the ground-state wavefunction of V^~\x) with energy eigenvalue 
= - 3{2A - 1) = (^ - 2)2 (case Hi of SUSYQM). 

We conclude that for Scarf I potential, the quadratic case leads to rather similar results to 
those obtained for the radial oscillator, the form of the rationally-extended potential becoming 
sensitive to the type of reparametrization made for the starting conventional one. 



3.2 Determination of wavefunctions 

The procedure used here to determine the wavefunctions of the rationally-extended Scarf I 
potentials being the same as that introduced in Section [2.21 we will only state the results. 



3.2.1 Linear case 

On choosing for V^'^\x) the Scarf I potential V^.B+iix), for instance, we obtain for the wave- 
functions of V^~\x) 

y/Eu p + a—[p — a)z 

where AAi+^ = AAi^'^~^\ = [y + A- B + \) {u + A + B - \), a = A- B-\, P = A + B-\, 
z = sinx, and 

= [;3 + « _ (^ _ c,)z] (^(1 + z)^ + /? + 1^ + (/3 - a)(l + z). 

The action of this first-order differential operator on the Jacobi polynomial Pi" ^'^~^^\z) can 
be easily proved to be given by 

0{«,/3)p(a-l,/34-l)(,) ^ 2(/3 -«)(/? + v)Plr^^\z), 

where P'^^\z) is a (z^ -|- l)th-degree polynomial defined by 

+ (/3 + a + 2.)-' (^p(-«)(,) - P<^f (.)) . (3.11) 

The partner potential wavefunctions (I3.10p can therefore be written in terms of -Pj"'f ^ {z) as 



, _^(l-sinx)5{^-^)(l+sinx)3{^+^) -(A-B-i.A+B-i), . , 

ihl >(x)=N,) '- ^ P,):, ^ ^Msmx), (3.12) 

^ 2A-l-2Bsmx v \ j 
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where 

(_) _ B f {2A + 2v)u\T{2A + v) C 

" ~2^\ {A-B+v+\) {A+B+v+\) r {A-B+v-\) T {A+B+v-\) ] ' ^^'^^^ 

Had we started from the Scarf I potential Va,b-i{x), we would have been led to computing 
the action of 

a^^'/^) = [/3 + a - (/? - a)z] ({I - z) A _ (a + 1)^ + (/3 _ «)(! _ z) 

on pj""*"^'^ Since Q^'^^ is changed into —Of'^^ under the permutation of a with 

combined with the transformation z — > — z, we would have obtained the same results (j3.12p 
and ()3.13p as before, up to some irrelevant overall sign. 

pIjj^i\z), V = 0, 1,2, . . ., form the second complete orthogonal set of polynomials with re- 
spect to some positive-definite measure that was constructed in [281 [29] by starting with some 
linear polynomial. As recalled in Appendix [Bl these Xi-Jacobi polynomials can be expressed 
as linear combinations of three classical Jacobi ones with constant coefficients. They are nor- 
malized in such a way that their highest-degree term is _2~'^~i(^'^"^°"'"^)2;^+^ (as compared with 
2-i'(2i/+a+/3^^!/ £qj. pj"'^)^^)^. In the same appendix, two special cases of well-behaved poten- 
tials V^~\x) outside the allowed range 0<-B<^ — lof parameter values are considered in 
connection with two so far unknown limiting properties of 

Observe that the first occurrence of the rationally-extended Scarf I potential (j3.5p and of the 
exceptional Xi- Jacobi polynomials in quantum mechanics can be traced back to [34], where 
it was also demonstrated that such a potential is shape invariant with a partner given by 

VA+l,B,cxt{x). 

We now plan to generalize Pj^'^P {z) to the more sophisticated extended potentials introduced 
in Section [3.1.21 

3.2.2 Quadratic case 

Let us start by considering the rationally-extended potential corresponding to case I and defined 
in equations ()3.6p and ()3.7p . In calculating its wavefunctions, we arrive at the first-order 
differential operator 



OT''> = V{z) 1^(1 + z)_ + /? + Ij - (1 + z)V{zl (3.14) 

where a = A — i? — ^,/3 = ^-|-i? — |, and T>[z) amounts to the function D[x) re-expressed in 
terms of z = sinx, while T){z) denotes its derivative with respect to z, i.e., 

V{z) = {(3- a- 2)[{f3- a- l)(/3- a- 2)z^-2{l3- a- l){f3+ a)z + {(3+ af+ (3- a- 2], 
V{z) = 2{P-a- l)(/3 - a - 2)[{f3 - a - 2)z - {(3 + a)]. 

Such an operator leads to a new family of (z^ -|- 2)th-degree Jacobi- type polynomials -Pi [/.^2('^)' 
defined by 

^(a,/3)p(„_i,/3+i)(^) = 4(/3 - a - l)(/3 - « - 2)2(z. + /3 - l)Pi'^%{z) 
= i^l""^ |to-«-l)(/?-«-2)(2z/ + /? + a)(^ + /3-l)z^ 
- 2{p - a - + a)[{i^ + P){2l' + /3 + a) + (3 - a - 2]z 
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+ + P + l){2u + p + a)[{p + af+(3-a-2] + 2{P-a- l)(/3 + af] P^"^^\z) 
+ A{l3-a-l){v + 0)[{[i -a-2)z-{(3 + a)]P^!f (z)}, (3.15) 

and the result for ipi ^ (x) reads 



^iKx)=My^ Plu+2 "{smx), 



where 

\2 



AA(-) = ^''-^1}^^-^\ {2A + 2u)u\ n2A + .)] 



X 



[{A- B + v + 1) [A + B + v + \) [A + B + v - \)] 



lM-l/2 



X [T {A - B + u - \)T {A + B + u - I)] 



-1/2 



For the case II potential given in equations (j3.6p and (jS.Sp . it turns out that the first-order 
differential operator 0^'^\ appearing in the calculation of its wavefunctions, satisfies the same 
type of symmetry relation with Q^'^^ of equation (j3.14p as that connecting O^'^^ with Q^'^^ 
in the linear case. In other words, O^'^^ can be obtained from Q^'^^ by simultaneously per- 
muting a with P and changing z into —z. As a result, no new family of Jacobi-type polynomials 
arises for such a potential, its wavefunctions being expressed as 

, ^,^.^(_)(l-sinx)^(^-S)(l+sinx)5{^+^) (A+B-i.A-B-i), . , 
iPl >{x) = {-iyMy^ Pi,u+2 ^\-smx), 

x[{A-B + u + ^){A-B + u-l) {A + B + i^+'l)y^^^ 
x[T{A-B + u-l)T{A + B + u-^)] , 

in terms of the {v + 2)th-degree Jacobi-type polynomials defined in p.lSp . 

Finally, for the case III potential defined in p.6p and (|3.9p . the counterpart of (|3.14p reads 

of ^ V{z) (^{l-z^)^^+P-a-{P + a + 2)z^ - (1 - z^)'D{z), 

where a, (3 are the same as before and 

V{z) = {l5 + a-l)[{(3 + a-2)z-{(3- a)f + {2a - 2){2(3 - 2), 
V{z) =2{(3 + a- l)(/3 + a - 2)[{f3 + a - 2)z - - a)]. 

This gives rise to a new family of Jacobi-type polynomials P^'^J^^{z), defined by 

= 8(/3 + a - 1)(/? + a - 2)2(/? + a + ,.)P^'^%{z) 

-{ {p + a + u + 2){p + a + 2u + l) 



1 



{f3 + a + u + l)(/3 + a + u + 2){P + a + 2u) 
X [{P - a){P + a) - {P + a + 2v){P + a + 2v + 2)z\V{z) 
+ [-{P + a + V + 1){P + a + 2v){P + a + 2v + 2) 
+ {P- afv + 2(/3 - a)v{P + a + 2u + l)z 
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+ {P + a + 2u){P + a + 2u + l){P + a + 2u + 2)z'^]V{z)\ P^'^'f^^z) 
+ 2{a + u){p + u)^{P + a + u + 2){P + a + 2u + 2)V{z) 
-[{f3-a) + {f3 + a + 2u + 2)z]'D{z)] P^fiz)}. 

The extended-potential excited-state wavefunctions can be written as 

(_) ^,(^) (l-sinx)^(^-^)(l + sinx)i(-^+^) (A-i?-l.A+i?-l) ^ . 
ru+i('')=K+i ^^^^ ^sV+s ''(smx), 

(-) _ {A-l){2A-3f ( {2A + U- 1){2A + 2u + 2)v\ T{2A + u + 2) \ ^''^ 
A/.+i - + ^)Y {A - B + v + l)T {A + B + v + 1) J ' 

It is obvious that the two families of Jacobi-type polynomials that we have just introduced 
are of a different nature. From physical considerations, the first set of polynomials Pi'^J%{z), 
whose lowest-degree one is quadratic in z (see Appendix [C]), is a good candidate for the still 
unknown complete, orthogonal family of X2-Jacobi polynomials. By analogy with classical and 
Xi-Jacobi polynomials, ^i"+2(-^) been normalized in such a way that its highest-degree term 
is 2~'^~'^ (^^'^'^^'^^^ z'^^'^ . By contrast, the other family P^°'j^^{z), starting with a cubic polynomial, 
cannot be complete. Note that the highest-degree term is now given by _2~'^~3 ^2!/+Q+/3+2^^;y+3^ 

Let us emphasize that for Jacobi-type polynomials, no splitting similar to that observed when 
going from lI^^-^{z) to L^^l_^_2{z) and l''^1_^_2{z) has been encountered. 

Finally, in cases I and II, the polynomial ^^^"2 ^' ^ (sinx) (or its counterpart) can be 

inserted in the extended potential ground-state wavefunction ipQ \x) to prove that the corre- 
sponding potential Va,b, cxt{x) is shape invariant, its partner being given by VA+i,B,cxtix) (see 
Appendix [D]). This again generalizes a result of [M] to the quadratic case. 



4 Final comments 

In the present paper, we have generated new exactly solvable rationally-extended radial oscillator 
and Scarf I potentials and we have constructed their bound-state wavefunctions. This has been 
made possible by generalizing a constructive SUSYQM method recently proposed in [27] and 
based on some reparametrization of the conventional superpotential, together with the addition 
of a rational term expressed in terms of a polynomial g{z), where z is some appropriately chosen 
function of x. The cases of linear and quadratic polynomials have been considered here. 

In the linear case, there appears a single rationally-extended potential, but it can be obtained 
by starting from an isospectral conventional potential with two distinct sets of reparametrized 
couplings. In contrast, the quadratic case leads to a variety of rationally-extended potentials, 
each of them being the partner of a single reparametrized conventional potential. Some potential 
pairs turn out to be isospectral, while in others, the rational potential has an extra bound state 
below the conventional potential spectrum. 

All rationally-extended potentials belonging to isospectral pairs have been demonstrated to 
be shape invariant as their conventional counterparts. 

Furthermore, by applying our SUSYQM approach, we have explicitly shown that the (z^-l-l)th- 
degree polynomials {v = 0, 1, 2, . . .) occurring in the bound-state wavefunctions of the extended 
potentials corresponding to the linear case are the Xi-Laguerre or Xi-Jacobi polynomials that 
were recently proved to form two sets of complete, orthogonal exceptional polynomials. We 
have then proposed several extensions of these polynomials valid for the quadratic case. Among 
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them, we have identified two different kinds of {u + 2)th-degree Laguerre-type polynomials 
and a single one of {u + 2)th-degree Jacobi-type polynomials, which are candidates for the still 
unknown X2-Laguerre and X2-Jacobi exceptional polynomials, respectively. 

Two interesting properties, dealt with in [27], have not been considered here, but are worth 
mentioning. To begin with, our first-order SUSY transformation may be combined with another 
one relating conventional potentials with different parameters to produce a reducible second- 
order SUSY transformation connecting conventional and extended potentials with the same 
parameters. In the linear case, this results in a second-order transformation admitting two 
distinct factorizations. 

The other point has to do with the possibility of discarding the restriction to real potentials, 
which has been implicitly made here. Considering also T^T-symmetric complex potentials indeed 
facilitates reconciling our approach to the requirement that the rationally-extended potentials 
be singularity free, hence generating so far unknown complex potentials with a real spectrum. 
The generalization of the present work along these lines is therefore an interesting topic for 
future investigation. 

When the present work was in its final stage, there appeared two interesting preprints, 
whose subjects are related to those considered here. In the first one [36], the existence of 
distinct factorizations of second-order SUSY transformations into products of two first-order 
ones, observed for the rationally-extended generalized Poschl-Teller potential of [27] and which 
would also be obtained here in the linear case, is discussed in the framework of type A 2-fold 
SUSY [22]. Necessary and sufficient conditions for such a situation to occur are derived and some 
relations to second-order parasupersymmetry and generalized 2-fold superalgebras are noted. 

In the second work [37], three infinite families of shape- invariant, rationally-extended radial 
oscillator, trigonometric and hyperbolic Poschl-Teller potentials are presented. They are ob- 
tained by deforming the corresponding conventional potentials in terms of their degree £ polyno- 
mial eigenfunctions and their bound-state wavefunctions are expressed in terms of Laguerre-type 
or Jacobi-type polynomials. In the radial oscillator and trigonometric Poschl-Teller cases, the 
first member of the infinite family, corresponding to ^ = 1, is shown to coincide with one of the 
potentials introduced in [3l] and re-obtained here in the linear case for the radial oscillator and 
Scarf I, respectively. 

The general expressions provided both for the potentials and the polynomials in [37] enable 
us to pursue the comparison with the results of [27] and those derived here. To start with, it 
can be easily seen that the first member (with i = 1) of the third family (hyperbolic Poschl- 
Teller) actually corresponds via some changes of variable and of parameters (x — > x/2, g — > 
B — A — 1, h — > B + A + 1) to the extended generalized Poschl-Teller potential constructed 
in [27] . Furthermore, the results corresponding to the second member (with i = 2) of the first two 
families agree with those of the present paper associated with the quadratic case and referred to 
as case I extended radial oscillator and case II extended Scarf I, respectively. It should be stressed 
that no case II extended radial oscillator (with its corresponding Laguerre-type polynomials) is 
found there. Whether the existence of such an alternative branch of potentials and polynomials, 
demonstrated in the quadratic case in the present paper, could be generalized to higher-degree 
polynomials g{z) would be an interesting topic for future investigation. 

A Examples of Laguerre-type polynomials 

In this appendix, we list the first few -£i°j_|_2i ^^^'^ ^s'u+s Laguerre-type polynomials and, 

for comparison's sake, also the corresponding classical and Xi-Laguerre polynomials. 

Laguerre polynomials 

4")(z) = i, 
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L^^>{z) = -z + a + l, 

L^^\z) = \[z'- 2(q + 2)z + (a + 2)(a + 1)]. 

Xi-Laguerre polynomials 

L^^\z) = -z-a-l, 
L^^\z) = z'^ -a{a + 2), 

L'f\z) = + (a + 3)2^ + a{a + 3)z - a{a + l)(a + 3)]. 

New Laguerre-type polynomials 

L'l^iz) =z^ + 2{a + 2)z + {a + 2){a + 1), 
L^^l{z) = -z^ - (a + 3)z^ + a{a + 3)z + a{a + l)(a + 3), 
L^"liz) = \\z^- 2{a + 1)(q + 4)^2 + a{a + \f{a + 4)], 
L^°^(z) = _^ 2az + a(a + 1), 

L^Q^iyz) = -z^ - (a - \)z^ + (a + 2)(a - \)z + (a + 2)(a + l)(a - 1), 
L^(z) = ^[z^ - 4z^ - 2{a + 3)(q - \)z^ + (a + 3)(q + 2)a(Q - 1)], 

(-2) = + 3(a - 1)-^^ - 3a(a - l)z + (a + l)a(a - 1)], 

l^^l^z) = \\z^ - ^az^ + 2{a - l)(3a + 4)z^ - 4(q + 2)a(Q - l)z 
+ (a + 2)(a + l)a{a - 1)], 
= jq[-z^ + 5(a + l)z^ - 10(a2 + 2a - l)z^ + 10(a + 3)(a + l)(a - l)z'^ 
- 5(a + 3)(a + 2)a(a - 1)2; + (a + 3)(a + 2)(a + l)Q(a - 1)]. 

B Limiting cases of extended Scarf I potentials 
and of Xi-Jacobi polynomials 

The purpose of this appendix is to review two cases where although the parameter B in the 
rationally-extended Scarf I potential p.Sp takes a value outside the allowed range < -B < ^ — 1, 
the potential remains physically acceptable and reduces to some known conventional potential. 
As a result, there exists a limiting relation between the wavefunctions of the former, given 
in equations (j3.12p and (j3.13p . and those of the latter, expressed in terms of some classical 
polynomials. The corresponding properties of the Xi-Jacobi polynomials will be demonstrated 
by starting from their known ones, proved in [29] • 



If we set i? = in equation ()3.5p , the sum of the two rational terms vanishes and we get the 
B ^ limit of the Scarf I potential, which is the one-parameter trigonometric Poschl-Teller 
potential Va,o{x) = A{A — l)sec^x. Its energy spectrum is given by equation ()3.ip and the 
corresponding wavefunctions can be expressed in terms of Gegenbauer polynomials as [38] 



iPi'^'^^ (x) = ATj^) (cos x)^Ci^^ (sin x) 



where 



( 



r{A)r{2A)ul{A + u) 



) 



1/2 



^T{A + ^)T{2A + iy) 
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Comparison with equations (|3.12p and (j3.13p leads to the relation 

hm (B - a) P("'^) (z) - (a + ^ + l)r(2a + l)r(a + u) A^+h). . 

The direct proof of equation (IB.lj) is based on the expansion of Xi-Jacobi polynomials as 
linear combinations of three classical Jacobi ones, 



13 - a{[i + a + 2v){l3 + a + 2v + 2) ^ 



{f3 + a + 2u){(3 + a + 2v + l) 



On letting /? go to a in {j3 — a)Pj^"'{'^(z), it is clear that the only surviving term is proportional 

to Py'^'^\z), which is known to be expressible in terms of Cy ^^•'(z) [32], thus leading to 
equation ()B.ip . 

On the other hand, we observe that although the Scarf I potential Va,b{x) is not defined for 
B = A—^ or A = 3+}^^ the same is not true for its extension VA,B,e^t{x), given in equation (|3.5p . 
The latter is indeed equivalent to the well-behaved conventional Scarf I potential Vj^_^^ A-^i^) 

with energy spectrum E^^'^^'^ = [A + v + 1)^, = 0, 1, 2, . . ., and wavefunctions ij^l ' ^''(x), 
obtainable from equation ()3.2p . For the same parameter values, we get from equation (|3.12p 



lim Vo \x) (X {1 — sinx) 4(1+ sin x)^ 4 [2A + 1 — (2^ — 1) sin x], 
B^A-\ 



which does not vanish for x ^ ^, hence is not physically acceptable. This explains the absence 

of an eigenvalue in the energy spectrum. The presence of the remaining eigenvalues, corre- 

^ ( i_ j^j^Q l.^ 

sponding to u = 1, 2, . . ., in ()3.ip hints at a limiting relation between Pj)_^i ^' ^ (-2) and 
Pyli {z) when i? ^ A — 2- In terms of a and /3, such a relation can be written as 



lim ef (.) = -^^(1 - ^fPt'^i^l . = 1, 2, . . . . (B.2) 



To prove equation ()B.2p . let us start from the a — > limit of the defining relation ()3.1ip of 
Xi-Jacobi polynomials. 



...... ,^|- V~ I ^u-l\^) 



r,(",/3) 



On using equations (22.7.18) and (22.7.15) of [3^, it is straightforward to transform the latter 
into 



limP^^f -_(/^±il±lLn _ .\ ^„p(l,/3)^.^_^,^^^p(l'/3)^ 



Another application of equation (22.7.15) then leads to the searched for result ()B.2p . 
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C Examples of Jacobi-type polynomials 

In this appendix, we list the first few Pi"j^2 ^^'^ ^Su+s Jacobi-type polynomials and, for com- 
parison's sake, also the corresponding classical and Xi-Jacobi polynomials. 

Jacobi polynomials 

P^^)(z),= i[(/3 + a + 2)z-(/3-a)]. 
Xi-Jacobi polynomials 

Pt'^\z) = ^(p^^i-if^ -a)z + (3 + a + 2], 

Pt^\z) = ^p\^^ {-if^ - + a + 2)^^ + P - a)^ + {(3 + a){(3 + a + 4)]z 

- iP-a)iP + a + 2)}. 
New Jacobi-type polynomials 

piT(') = jtr ^[(^ - " - - « - 2)^' 

4(p — a — l)(p — a — 2) 

-2{(3-a-l){[i + a + 2)z + {(3 + a + 2f + [i-a-2], 

^a-"<^> = 8(ff-„-iKff-„-2) «^ - " - - " - ^'c' + ° + 

- (/? - a - - a)(/3 - a - 2) + 2(/? + a)(/3 + q + 4)]^^ 

+ (/? + a + 2)[(/3 - a - 2)(2/? - 2a + 3) + (/3 + a)(/3 + Q + 4)]z 

- (/? - «)(/? + a + 2)2 - (/? - a - 4)(/3 - Q + 2)}, 

+ 3(/3 + a)(/3 + a - l)(/3 - a)z^ - 3(/3 + a)[(/3 - a)^ + /? + a - 2]z 

\2 



+ iP- a)[{(5 - ay + 3(3 + 3a- 4]}, 
1 



16(/3 + a - + a - 2){p + a + 1) 
X + a - 2)(/3 + a - l)(/3 + a + + a + 4)z^ 
+ 4(/? - a)(/3 + a - + a + l)(/3 + a + 2)z^ 
-2{p + a + l)[(/3 - a)2(3/3 + 3a + 2) + (/3 + a - 2)(/? + a + 4)]^^ 
+ 4(/3 - a)(/5 + a + !)[(/? - af + (3 + a - 2]z 
- (/? - a)^ - 2(/3 - a)2(/? + « _ 4) + (/? + a - 2)(/? + a + 4)}. 

D Shape invar iance of rationally- extended potentials 

In this appendix, we prove the shape invariance of the rationally-extended radial oscillator 
potentials, defined in equations (12.150 and ()2.17p (the latter with upper signs only), as well as 
that of the rationally-extended Scarf I potentials given in equations ()3.6|) . (|3.7p . and (|3.8p . 
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In such a picture corresponding to case i of SUSYQM, the rationally-extended poten- 
tial is considered as the starting potential V^~^\x) and its partner is determined from equa- 
tion (j2.6p as V^~\x) = V^^\x) + 2W'{x), where the superpotential is now given by W{x) = 
—d[liii^^\x)]/dx, with 'ijjl^\x) = ipQ 

On using the expressions found for ^ (x) in Section 12.2.21 and the results of Appendix |X1 
it is straightforward to show that for case I and II extended radial oscillator potentials, 

W(x) = Wi(x) + W2(x), Wi(x) = -u;x-^-^, 

2 X 

and 

= i (u;x^ + 2/ + 3)^-2(2/ + 3) " {ux- + 21 + 5)^ - 2(2/ + 5) J ^^"^^ 

~ / uox'^ + 21-l Lox'^ + 21 + l \ 

W2{x) = Aux \ , TT - f..^-?. , o, , , o.o, , ^^ (case II), 



(u;a;2 + 21- 1)2 + 2(2/ - 1) {ujx"^ + 21 + 1)2 + 2(2/ + 1 

from which it directly follows that 

2W'{x) = -Vi^eA^) + V^+i,ext(^) + u:. 

Hence, in such cases, the partner of Vi^cyA.{x) is V/+i^oxt(3^) +'^) which proves the shape invariance 
of the former. 

Similarly, from Section 13.2.21 and Appendix [Cl we obtain for case I and II extended Scarf I 
potentials, 

W{x) = Wi{x) + W2{x), Wi{x) = ^tanx - ^secx, 

and 



W2{x) = 2{2B - 1){2B - 2)cosx 

' (2B - 2) sinx - {2A - 1) {2B - 2) sin x-{2A + 1) 



Da,b{x) Da+i,b{x) 
W2{x) = 2(2S + l)(2B + 2)cosx 

■ {2B + 2) s\nx - {2 A - 1) {2B + 2) sinx - {2A + 1) 

Da,b{x) Da+i,b{x) 



(case I), 



(case II), 



where Da,b{x) denotes the denominator function D{x) of equations (|3.7p and (j3.8p . associated 
with some specified parameters A, B. Since 

2W'{x) = -VA,B,eAx) + VA+l,B,eAx), 

the partner of VA,B,cxi{x) is VA+i,B,cxt{x)^ thus completing the proof. 
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